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We explore the energy dynamics of a two dimensional (2D) trapped Bose-Einstein condensate 
(BEC) excited by a moving red-detuned laser potential (RDLP). The trap is a combination of 
a power-law (PL) potential cutoff by a hard- wall box potential (HWBP). It is found that by a 
restricted measurment of the energy within the boundaries of the HWBP, the energy demonstrates 
oscillations indicative of solitons. It is then demonstrated, that the geometry of the PL potential 
influences the lifetime of these oscillations, i.e., the lifetime of the ensuing solitons inside the HWBP. 
We argue that the energy dynamics in this setup are a good tool for measuring their lifetime. It is 
also found, that the condensate density dynamics display oscillatory patterns of a magnitude and 
order controlled by the velocity of the RDLP. A connection to oscillations in the chemical potential 
dynamics is also discussed. Essentially, we suggest future experiments for this kind of setup, which 
would measure the phonon energy dynamics to trace the lifetime of solitons. 

PACS numbers: 03.75.-b,67.85.De,67.85.Jk 



I. INTRODUCTION 

The investigation of a Bose-Einstein condensate 
(BEC) excited by an obstacle has received considerable 
theoretical and experimental interest in the last few years 
[THIS]. So far, mostly repulsive obstacles have been cre- 
ated experimentally by a blue-detuned laser beam [5lU6| . 
Recently, however, the use of attractive obstacles created 
by a red detuned laser beam has become widespread, ex- 
perimentally [TTH2"5] as well as theoretically |2ftH31| . This 
red detuned laser beam creates a negative potential well, 
also known theoretically [TTJ [27J EI] an d experimentally 
OH E3J E2] as a "quantum tweezer" (QT) or "dimple- 
trap" (DT). DTs have become important tools in quan- 
tum computation [31 j . The importance of understand- 
ing the dynamics of BEC trapped in a moving DT, its 
response to the velocity of the DT and its depth, lies in 
future applications of atom qubit transport in an optical 
lattice (OL) to implement collisional quantum computing 
gates [31] . They have also been used to transport a BEC 
for long distances [23] , for BEC production [211 [24] , and 
the control of the scattering length [T7j. Earlier, Carpen- 
tier et al. |17j investigated the extraction of atoms from a 
cigar-shaped BEC using a laser tweezer. Diener et al. [27] 
proposed a QT for extracting a number of particles from 
a harmonic trap. This QT was a moving dimple, and the 
extraction of atoms was by Landau-Zener tunneling. 

In this paper, we consider a 2D hard-sphere Bose gas 
confined by a general power-law trap, cutoff by a hard- 
wall box potential (HWBP) boundary. It is the same 
trapping geometry we used earlier [33], except that we 
add here to the external trap a moving DT potential 
instead of a repulsive obstacle. 

The chief goal of this paper is to explore the effects 
of trapping geometry and dimple speed on the soliton 
lifetime and condensate density dynamics, the most im- 



portant finding being the observation of oscillations in 
the system's energy which can be used to trace the life- 
time of the solitons. Previously Parker et al. [34] demon- 
strated energy oscillations of solitons as well, by restrict- 
ing their calculations to the dimensions of the dimple 
trap. Similarly, we restrict our calculations to the area 
of the 2D HWBP in order to capture the dynamics of soli- 
tons. Moreoever, oscillations in the condensate density 
dynamics are also observed. 

In our previous publication [33j . it has been demon- 
strated that a blue-detuned laser potential (BDLP) mov- 
ing inside a 2D power-law (PL) trap plus HWBP bound- 
aries can lead to fascinating self-interfering matter wave 
(SIMW) patterns. In this work, we invert the "sign" 
of the BDLP to make it attractive and investigate the 
ensuing BEC dynamics thereof. As demonstrated be- 
low, the moving DT traps a BEC fragment and splits it 
away from the central BEC, transporting it towards the 
HWBP boundary. When the dimple reaches the HWBP, 
the soliton is reflected from the hard wall as the DT 
leaves the trap. Once the soliton bounces of the hard 
wall, it heads back towards the central BEC and merges 
again with it. In fact, the splitting and merging of BECs 
has been performed earlier to observe many fascinating 
phenomena [hH I35H35] . For example, Scherer et al. [T3] 
successfully observed vortex formation by splitting a har- 
monically trapped BEC into three sections using a blue- 
detuned laser, and then allowing those sections to merge. 
Chang et al. [38] has demonstrated experimentally, us- 
ing a BDLP, that solitons and dispersive shock waves 
(DSWs) can be produced by merging and splitting of 
BECs, respectively. Here in this work, dispersive waves 
(DWs) are also produced by merging BECs, but using a 
red-detuned laser potential (RDLP) or DT. 

Our key results are (i) the lifetime of solitons gener- 
ated by the above merger of two BECs depends on the 
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shape of the PL potential; (ii) the condensate density 
oscillates with time due to oscillations in the chemical 
potential, and the oscillatory pattern depends greatly on 
the velocity v of the dimple as it becomes more and more 
ordered approaching sinusoidal behavour with increasing 
v; (iii) for very small v approaching zero, the RDLP be- 
haves like a BDLP, as it can generate the same cross-like, 
SIMW patterns similarly to the ones observed earlier [33J • 
To the best of our knowledge, our current findings per- 
taining to the energy dynamics in a trapped BEC excited 
by a QT have not been reported earlier. 

Our paper is organized as follows. In Sec(ll] we give 
a brief description of the method; for more details and 
information please refer to our previous article |33| . In 
Sec lIIII we demonstrate our results with focus on the en- 
ergy oscillations, and in Sec|IV|we summarize our results. 



II. METHOD 

We consider again N hard core (HC) bosons inside a 
2D PL trap cutoff by a HWBP boundary [33J . The model 
we use is the same as in our previous work [33]; except 
that we apply here a moving RDLP instead of a BDLP. 
It is the same as the BDLP, but with a negative ampli- 
tude. Consequently, it is a trough instead of a barrier. 
The RDLP continues its motion and leaves the HWBP 
boundary without returning to it. 

The split-step Crank Nicolson method [39] is applied 
to solve the 2D time-dependent Gross-Pitaevskii equa- 
tion (TDGPE) [Eq.(l) in Ref.[33]]: 



presented here are in the transient stage of the simula- 
tion, i.e., after and not including the initialization pro- 
cess. (Refer to [33 for more specific details on the ini- 
tialization process.) The energy E(t) (function of time) 
is computed using Eq.(7) in Ref . [33] 



E(t) 
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\Vip(x,y;t)\ 2 + V(x,y;t)\cp(x,y;t)\ 2 + 

(3) 



where the limits of the integration are only over the area 
of the HWBP, i.e., -10 < x < 10 and -10 < y < 10. As 
a result of this limitation, we largely capture the energy 
dynamics of the solitons inside the HWBP, and E{t) will 
not be constant with time. As the DT begins to move, 
starting from the center of the trap, it pumps energy into 
the system as it "climbs up" the PL trap. Once it leaves 
the HWBP, it is not anymore included in the integral 
(3). The N , A, v, /3, E(t), and time t, all have the same 
units as before [33] : The lengths and energies are in units 
of the trap = \Jfij (2mwf, ) and hujh , respectively. 
A is in units of huiho, P in (a?io)~ 2 , v in ai lo , t = rw/ lo 
is unitless, N is in units of (v^a 2 ^) -1 , and the density 
\(p(x,y;t)\ 2 is in (a/j ) -2 . The anisotropy parameter k is 
unitless. It must be further noted, that the wavefunction 
tp(x,y;t), as well as its gradient Vtp{x,y;t), are assumed 
to be zero at the HWBP boundary. This is in order to 
enforce the hard wall effect which is of infinite potential 
height. 
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where, as before, V(x,y;t) — V(x,y;t)/huiho is an ex- 
ternal potential given by the combination 



V{x,y;t) = 

Jd^r +*\yn + AeM-Pl* 2 + (y~vt) 2 )]}, 

(2) 

Who being the trapping frequency, Af the coupling con- 
stant, a the strength of the 2D PL trap, pi and p2 any 
positive numbers, k — cu y /uj x the anisotropy parame- 
ter, A < the amplitude of the RDLP, /3 the exponent 
determining the width of the RDLP, and v its veloc- 
ity. (p(x, y; t) is the wavefunction of the system, where 
J^dxJ^dy\^(x,y;t)\ 2 = 1. As in Ref.[33], the sys- 
tern is initialized by a gradual introduction of the nonlin- 
earity and a moving RDLP whose amplitude and velocity 
are gradually increased with time until they reach their 
desired values in this simulation. However, the results 



III. RESULTS AND DISCUSSION 

A. Splitting and Merging of the BEC 

The atoms are attracted to the well of the RDLP, 
thereby creating a density peak as the atoms fill the 
trough. In other words, the DT traps a BEC fragment 
which is extracted from the main part of the reservoir. 
As the DT moves, it splits and carries the fragment away 
from the main part towards the HWBP boundary. When 
the DT fragment reaches the hard wall, it bounces off 
and returns in the opposite direction as the RDLP leaves 
the HWBP area. This is demonstrated in Figs.[T and 
[2] which display the 2D BEC densities \ip(x,y;t) 2 at 
different times t = 4.5 and 5.4, respectively. The sys- 
tem is a trapped BEC in an isotropic 2D harmonic trap 
(pi=P2=2) cut off by a HWBP boundary, as in Ref. [53] . 
The hard walls are located at x, y—± 10, as displayed 
in Figs.[l] and [2] The parameters used are JV = 10, 
A = -30, P = 4, v = 2, and k = 1. In Fig. [l] the 
separated BEC fragment transported by the DT collides 
with the hard wall at y = 10 and t — 4.5. In Fig. [2] 
the BEC fragment has bounced off the wall and is on 
its way back to the center of the trap, when the DT has 
already left the HWBP. Since the BEC fragment is no 
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M = 10, A = -30, t = 4.5 

/3 = 4, v = 2, k= 1, pi =p 2 = 2 
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FIG. 2: (Color online) As in Fig.[T] but for t = 5.4. The units 
FIG. 1: (Color online) Two dimensional (2D) density are as in Fig.[l] 
\<p(x, y\ t)\ 2 at t = 4.5. The system is a hard core Bose gas 
in a 2D harmonic trap plus box-potential boundaries [33] ex- 
cited by a red-detuned laser potential with N = 10, A — —30, 
/? = 4, v = 2, and k = 1. A is in units of huiho, ft in (ah )~ 2 , 
v in a^o, t = rcj^o is unitless, Af is in units of (v20h ) _1 , and 
y; t)\ 2 is in (ah ) -2 . x and j/ are in units of aho- 



v = 0.2. t = 1.0 



longer trapped by the DT, it expands in width and it 
declines in amplitude as it merges into the central part 
of the BEC. By monitoring the energy dynamics of this 
system in the transient stage of the simulation, one ob- 
serves an oscillatory pattern indicative of soliton excita- 
tions. According to Frantzeskakis et al. 40 and Weller 
et al. [41j . the merger and interference of two conden- 
sates leads to the generation of dark solitons. (See the 
upcoming Sec III C ) 

Previously, Dicner et al. [27] allowed the DT to travel 
far away from the trap center, until no atom was left 
in the dimple. In contrast, although our dimple trap 
keeps moving and eventually leaves the trap, it looses all 
the BEC at the HWBP boundary. Hence, an important 
feature in our work here is that we explore the effects of 
adding a HWBP to the quantum tweezer system. 



\<p(x,y;t)\ 




FIG. 3: (Color online) As in Fig.[T] but for 
t = 1.0. The units are as in Fig.[l] 



0.2 and 



B. Blue-laser like effect 

For very low velocities v < 1, the DT produces the 
same effects as a BDLP. In this regard, the mass trapped 
inside the moving DT generates circular dispersive waves 
(DWs) which expand and propagate away from the dim- 
ple. The outgoing DWs interfere with incoming DWs 
reflected from the hard wall. This is demonstrated in 
Fig|3]for the same system as of Fig. [l] but for v = 0.2. 
The system of Fig. |4j is for an even smaller velocity of 
v =0.01 at time t = 3.0, where in this case it displays a 



cross-like SIMW interference pattern resembling that of 
a Fraunhofer diffraction pattern. The later was observed 
earlier for a BDLP [33]. For further details on the BEC 
dynamics of Fig. [3j we refer the reader to the movie 
clip in the supplementary material [xxxx]. One could 
argue that the mass trapped inside the dimple behaves 
like a repulsive obstacle at these low velocities. This is 
because its mass and shape are almost stabilized by the 
attractive forces of the DT once it is separated from the 
central BEC. Therefore, an important finding here is that 
the RDLP with very small v behaves like a BDLP. 
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v — 0.4, P — 2 (open squares); and v — 2, /3 — 4 (open circles). 
The arrow indicates the onset time for solitonic nucleation, 



FIG. 4: Same as Fig. [3j but at t=3.0 for v=0.01. t ons . E(t) is in units of Tiloho and the rest of the units are as 

in Fig. [I] 



C. Energy dynamics 

In this section, we explore the energy dynamics of 
our systems. The goal is to explain their basic features 
before going into their implications. The basic feature 
is that the energy dynamics reveal an oscillatory behav- 
ior as demonstrated in the upcoming figures. Our chief 
argument, is that one can use this behavior to monitor 
solitons and their lifetime if one could measure E(t) ex- 
perimentally. We therefore suggest future experiments 
to measure the dynamics of the phonon intensity contin- 
uously emitted and reabsorbed by the solitons to trace 
their energy dynamics. 

In fact, Parker et al. [32] investigated earlier the dy- 
namics of a soliton oscillating in a tight dimple-trap em- 
bedded in a weaker harmonic trap. His system was quasi- 
one dimensional, and he found a characteristic beating 
period for his solitons. Parker et al. computed the en- 
ergy of his soliton by integrating the energy functional 
over the width of his dimple trap only. In contrast, our 
system is two-dimensional, and we investigate (3) over 
the area of a box enclosing solitons and phonons. Ear- 
lier, it was reported that soliton frequency is affected by 
dimensionality, the presence of other solitons, as well as 
trap anisotropy [4T1 14"51 134] . 

Figure [5] displays the oscillations in the energy as 
a function of time, E(t), for the same system of Fig. 
[I] but for different v and j3: (open triangles) v = 0.3, 
,3 = 2; (open squares) v — 0.4, (3 — 2; and (open 
circles) v = 2, (3 = 4. These oscillations are similar 
to those in Figs. 11-13 of our previous article [33] , where 
a BDLP was used. However, whereas in the case of a 
BDLP the energy begins to oscillate once the BDLP is 
moved with any velocity, in the case of an RDLP the 
energy oscillations suffer a delay for a certain time for 



some velocities such as for v = 2 [Fig. [5] (open circles)]. 

For v > 0.3, the initial energy of the system from t = 
to t ~ 2 in Fig.[5] is negative, because for the parameters 
chosen the negative RDLP energy outweighs the combi- 
nation of kinetic, bosonic interaction, and external PL 
trap energies. In contrast, the E(t) for v = 0.3 is posi- 
tive at all times as the system has such a kinetic energy, 
which in combination with the interaction and trap en- 
ergies outweighs the RDLP depth. It could be argued, 
that for v > 0.3, more material is initially trapped in the 
DT reaching negative potential energies. Consequently, 
the velocity of the dimple-trap manipulates the sign of 
the initial effective scattering length of the system. As 
time progresses, E(t) increases and becomes positive for 
v > 0.3 in Fig. [5] as the moving obstacle transfers en- 
ergy to the system in the initial stages of the simulation 
[45] , As the DT moves up the PL trap, the energy of the 
system increases since the energy levels inside the dim- 
ple increase with position |27j . This continues until the 
dimple leaves the system, at which point the reflected 
fragment merges with the main BEC part ("reservoir") 
and solitons are formed. Eventually, E{t) displays a tran- 
sient similarly to the one reported earlier by Law et al. 
[45] for a moving repulsive obstacle. However, whereas 
in Ref.j45] the system is a homogeneous BEC, where the 
energy stabilizes after the transient, in our case the tran- 
sient is followed by oscillations. Next to this, when the 
RDLP leaves the system, the effective local interaction, 
in the neighborhood of the RDLP, jumps from being at- 
tractive to repulsive. The RDLP hence changes the sign 
of the scattering length locally in the vicinity of dimple- 
trap |17) . By moving the laser beam, the region with a 
"local" negative scattering length moves with it. 

Thus, the principle feature in Fig. [5] above, is that a 
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FIG. 6: Energy dynamics for the same system as in Fig.[T] 
but for «=0.01 (open circles), and 0.20 (solid squares). 



transient is observed using an RDLP and that the initial 
scattering length is DT- velocity dependent. The oscil- 
lations in the energy following the end of the transient 
are indicative of dark solitons in the BEC which are a 
consequence of matter-wave interference of two merging 

BECs [SJIM1H]- 

It can be further argued, that before the merger the 
initial distance between the DT and the central part of 
the BEC is within the range where the interaction energy 
exceeds the kinetic energy. As a result, a dark soliton is 
nucleated due to nonlinear interference [H]. As the DT 
leaves the HWBP, the originally trapped BEC inside the 
dimple bounces off the hard-wall. Consequently, it starts 
to expand creating a shock wave which upon interfer- 
ence with the central BEC excites solitons. This leads us 
to the conclusion, that the rapid change in the effective 
interaction as the DT leaves the HWBP, is larger than 
the quantum pressure at the edges of the trap 50J . The 
shock wave pushes the main part of the BEC situated 
around the center of the PL trap slightly backwards. As 
such, BECs "communicate" together by this shock wave. 
The backflow around the released BEC fragment is still 
another way in which these two parts can communicate 
together [51] . The released BEC fragment temporarily 
acts like an (expanding) BDLP obstacle by itself, until it 
merges again with the central BEC part. 

For very small velocities, the energy dynamics in the 
initial stages of the RDLP motion behaves just like that 
for a BDLP. Figure [6] compares the energy dynamics for 
a very slow RDLP with that of a much faster one. The 
system is the same as that of Fig.|4| Open circles: v = 
0.01; solid circles: v = 0.2. For v = 0.01, there are 
initial chaotic oscillations similar to those in Fig. 13 (solid 
circles) of Ref. [33j . For v — 0.2, the oscillations follow an 
ordered pattern as they "climb up" the PL potential well. 
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FIG. 7: Energy dynamics for the same system as in Fig. [I] 
but different anisotropies k: 0.5 (crosses); 1.0 (times); 2.0 
(starred times); 4.0 (open squares); 6.0 (solid squares); 8.0 
(open circles); and 10.0 (open triangles). 



D. Soliton lifetime 

In this section, the effects of the trapping geometry on 
the lifetime of solitons is demonstrated. As noted earlier, 
the energy oscillations following the transient in Fig. [3] 
are indicative of dark solitons. It is demonstrated next, 
that for certain trapping geometries these oscillations can 
have a long life time; but for others they are damped 
within a certain period of time beyond the end of the 
transient, when the DT is about to leave the HWBP area. 
The decay in these oscillations signals that upon a merger 
of the two separated BECs, the solitons become unstable 
and decay in the high density buldge formed during the 
merger |38j . Particularly, it is found that the lifetime of 
the solitons can be extended by increasing the anisotropy 
k of the trap within a certain range, less than some k c , 
and that solitons vanish completely for k > k c . Here, k c 
can be considered some value beyond which solitons are 
suppressed, which is essentially system-dependent. 

Figures [7] and [8] display E{t) for the same system of 
Fig. [l] but various anisotropies k and PL-trap exponents 
Pi and p2- Increasing the anisotropy k between 0.5 and 
4.0 [FigJT^a)] yields a longer lifetime for the solitons. 
However for k > 4 [Fig{7]jb)], solitons are suppressed 
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FIG. 8: Energy dynamics for the same system as in Fig.[T] 
but different exponents pi = P2: 2.0 (crosses); 4.0 (times); 6.0 
(starred times); 7.0 (open squares); 8.0 (solid squares). 



FIG. 9: Condensate density dynamics for the same system 
as in Fig. IT] but different u=0.2 (open circles); 0.22 (solid 
circles); 30 (open squares); and 40 (solid squares). 



since there are no ordered oscillatory patterns anymore. 
Note that a noisy energy pattern is observed instead. The 
system in this case is in a chaotic state. 

On the other hand, tightening the trap by making 
Pi and P2 larger, only small solitonic lifetimes can be 
achieved between p\ — p2 € [2,6] after which solitons 
are suppressed. Yet again, the system displays a noisy 
energy pattern for pi — p 2 > 6. It has been indicated ear- 
lier [32], that solitons in extremely elongated geometries, 
are amenable to exhibit longer lifetimes since transverse 
excitations are frozen out [52] ■ However, this is under 
the condition that the transverse condensate size is of 
the order of the healing length. Basing on this, the sup- 
pression of solitons in Fig{7]jb) for k > 4 is because the 
transverse condensate size has become smaller than the 
healing length of the condensate. 

It should also be noted, that in a confining geom- 
etry such as our trap, cut off by a HWBP, the sound 
(phonons) emitted by the solitons cannot escape. Conse- 
quently, the noise observed in Figs(7^b) and [8] is largely 
due to phonons. Earlier, it has been demonstrated that 
dissipation of solitons is induced by deviations from har- 
monic trapping |42|, 153) . In addition, trap anisotropy 
breaks the soliton-sound equilibrium leading to a faster 
decay [33]. Indeed, our results are in line with the latter 
two findings. It has also been shown that, a single barrier 
in a Gaussian trap leads to soliton decay [54] . 

E. Condensate Dynamics 

In this section, it is shown that the BEC density dy- 
namics displays oscillatory behavior of a magnitude and 
order controlled by the velocity of the DT. In addition, 
the effect of trap anisotropy is demonstrated on the be- 
havior of Pft(0) 0;t). Oscillations in p FT (0,0;t) have 
been reported earlier [33] for the case of a moving BDLP. 
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FIG. 10: Chemical potential dynamics for the same system 
as in Fig.[l] but for different f=0.2 (open circles); 0.22 (solid 
cirbles); 30 (open squares); and 40 (solid squares). 



But in a DT, the bosons have the ability to condense fur- 
ther into the lowest energy level of the RDLP in contrast 
to the BDLP, where condensation is only into the lowest 
energy level of the PL trap 
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FIG. 11: (Color online) Condensate density dynamics 
PFr(0,0;i) for the same system as in Fig.[l] but for differ- 
ent values of k: 0.5 (crosses); 1.0 (times); 2.0 (starred times); 
4.0 (open squares); 6.0 (solid squares); 8.0 (open circles); 10.0 
(solid circles). 



Figure [9] displays the dynamics of pi?r(0,0;f) for the 
same system of Fig{T] but different DT velocities v. Open 
circles: v—Q.2; solid circles: v—0.22; open squares: u=30; 
and solid squares: w=40. As the velocity is increased, 
Pft(0-> 0;t) decreases in general, as there is an expected 
time-scale for condensate formation When v is in- 
creased the BEC has less time to grow causing a reduc- 
tion in its fraction. There seems to be a sudden rise in 
/9ft(0,0;£) for t ~ 10 and v = 0.22 (solid circles), in 
contrast to v = 0.2 (solid circles). For v of the order ~ 
0.1, the oscillations are not regular and chaotic, but for 
v of order ~ 10, the oscillations are quite regular and 
close to being sinusoidal. A reason for the observed BEC 
density oscillations can be given following Diener et al. 
[27] , who proposed a DT for extracting a desired number 
of atoms from a BEC in a harmonic trap. They noted, 
that while their moving DT is inside the BEC, the inter- 
action between the harmonic trap and the DT is stronger 
than the self-interaction energy of the bosons inside the 
DT. Consequently, the fluctuations of the number of par- 
ticles in the DT is large. The BEC density oscillations 
are therefore a result of a periodic process of particle ex- 
citations out of, and particle deexcitations back into the 
BEC of the DT. The bosons are therefore able to oscil- 
late between the condensate ground state and the excited 



FIG. 12: (Color online) Chemical potential dynamics fj,(t) for 
the same system as in Fig. [I] but for different values of k with 
the same labelling as in Fig. 1 1 1 [ 



states of the PL trap. 

Another reason for the observed oscillations can be 
stated as follows. The BEC inside the DT is formed when 
the chemical potential fi(t) of the system approaches the 
ground-state energy of the DT. When the DT is set into 
motion, it "climbs" up the PL potential barrier. Con- 
sequently, the ground-state energy of the DT increases 
with time. Therefore, when eg(t) (say) at some time ap- 
proaches fi(t), a BEC begins to grow. When fi(t) dis- 
tances itself from eo(i), the condensate begins to "evap- 
orate" back to the PL trap —the main reservoir, so to 
speak. 



Figure 10 displays the dynamics of fi(t) for the same 
systems as of Figj9] One can see that fj,(t) for v of the 
order of 0.1 is much larger than v of the order 10. Since 
fi(t) is oscillatory, one observes BEC density oscillations. 
Figures [TT] and [12] display the effect of the anisotropy on 
Pft(0, 0;£) and fi(t), respectively for the same system 
of Fig[l] and a range of K from 0.5 to 10. In general, 
Aft(0, 0;i) seems to decline if k is increased, but the 
oscillations remain. In contrast, /j,(t) rises with increasing 
k and its oscillations are suppressed. 

Earlier, Uncu et al. [53] investigated the properties of 
a BEC in a harmonic potential plus an eccentric static 
DT. It was found, that by increasing the distance between 
the DT and the center of the harmonic trap, the BEC 
fraction declined. In contrast, our work demonstrates 
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that the BEC fraction (which is proportional to the BEC 
density) oscillates under a continuous displacement of the 
DT. However, one must note that Uncu et al. used a 
Dirac-Delta function DT, whereas our DT is a Gaussian 
with some finite depth. It could therefore be argued, that 
a different shape for the DT potential can lead to different 
BEC dynamics. Next to this, in a tight stationary DT, 
suddenly introduced to an atom chip [53] , the condensate 
number of atoms was observed to increase with time. 
However, for our moving DT, the condensate number of 
atoms seems to oscillate instead. Thus, a moving DT 
does not allow the condensate to grow monotonically in 
a 2D trapped Bose gas as in our case. 

Finally, Stamper-Kurn et al. [20] explored the forma- 
tion of a BEC in a DT. They demonstrated repeated con- 
densation in the DT, i.e., oscillations in the condensate 
fraction with time similar to the oscillations observed in 
our systems. 

IV. CONCLUSIONS 

In summary, then, we have explored the energy dy- 
namics of a 2D BEC trapped inside a PL potential cut- 
off by a HWBP and excited by a moving DT. It was 
demonstrated that an attractive obstacle, the DT, leads 
to peculiar features in the energy dynamics, not observed 
in the case of a repulsive obstacle [33]. Particularly, for 



very small velocities of the DT, the RDLP behaves like 
a BDLP, causing again the generation of DWs leading 
to self-interfering matter waves. Further, the lifetime of 
solitons generated inside the HWBP by the method out- 
lined in Sec. Ill is largely influenced by the geometry 
of the PL trap, i.e. its confining strength (j>\ and P2) 
and anisotropy k. An increase in k up to a certain value, 
say k c , causes the lifetime of the energy oscillations to 
increase. But beyond a certain n c (which is system de- 
pendent), the oscillations are suppressed. On the other 
hand, an increase in pi and P2 leads only to a suppres- 
sion of the energy oscillations. Note that as one tightens 
the trap, the area of the BEC is reduced. Similarly, an 
increase in k, compresses the BEC in the direction of DT 
motion, away from the HWBP. Consequently, there is less 
interaction between BEC and HWBP. Next to this, the 
condensate density [pi?x(0, 0; t)] dynamics displays oscil- 
lations related to oscillations in the chemical potential 
n{t) of the system. For an isotropic trap, the oscillations 
in pft(0, 0; t) are chaotic for low v, but become more or- 
dered as v is increased. As one increases k, oscillations 
in pft(0, 0; t) remain, but disappear for p(t). We believe 
that our paper will stimulate experimental and theoreti- 
cal studies of trapped BECs inside a HWBP excited by a 
DT, particularly since there is a lack of literature in this 
regard. In the future, we shall explore the dependence of 
the critical parameter k c on the system parameters. 
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